#include "./mds_density_test.h"


void mds_density_box( const mds_frame *frame, mds_density *little_box )
{
	assert( little_box->ix > 0 );
	assert( little_box->iy > 0 );
	assert( little_box->iz > 0 );
	
	little_box->itaille = frame->box.x / 4;
	
	
}


void mds_density_number( const mds_frame *frame, mds_vector *r, mds_density *little_box )
{
	for(i=0; i<frame->num_particles; i++)
	{
		little_box->ix = floor( r->x / little_box->itaille );
		little_box->iy = floor( r->y / little_box->itaille );
		little_box->iz = floor( r->z / little_box->itaille );
		
		density [little_box->ix][little_box->iy][little_box->iz] += 1;
	}
	
	
	
	
	
}
